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Abstract. — Random matrix models of disordered bosons consist of matrices 
in the Lie algebra g = sp„(M). Assuming dynamical stability, their eigenvalues 
are required to be purely imaginary. Here a method is proposed for constructing 
ensembles {£,P) of G- invariant sets £ of such matrices with probability measures 
P. These arise as moment map direct images from phase spaces X which play an 
important role in complex geometry and representation theory. In the toy-model 
case of n = 1, where X is the complex bidisk and P is the direct image of the 
uniform measure, an explicit description of the spectral measure is given. 



The goal of this article is to point out a method for constructing ensemble prob- 
ability spaces in a bosonic setting. We do this in the context of the model in 
|LSZ] where due to the requirement of dynamical stability the generators 



of the time evolution are symplectic matrices of elliptic type, i.e., their eigenvalues 
come in pairs ±iAj with Xj > 0. We denote the set of these elements of the 
symplectic Lie algebra q = sp2„(M) by £. 

To put this in perspective it is perhaps useful to recall that in the case of fermions 
the basic model spaces consist of Hermitian operators H such that iH is, for ex- 
ample, an element of one of the Lie algebras su„, so„ or usp2n ol the classical 
compact groups. For example, in the case of su„ the basic model space becomes 
the vector space of all Hermitian matrices. In that case the model probability 
distribution can be chosen to be Gaussian, up to constants the density being 
f{H) := e-'l^"'. Since one is usually most interested in associated eigenvalue dis- 
tributions, it is entirely appropriate that the model densities are invariant with 
respect to the action by conjugation of the compact group K at hand, i.e., the 
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adjoint representation of K. In the case of sUn this is just the action of the spe- 
cial unitary group K = SU„. Reducing the symmetry by this action, one obtains 
associated classical Weyl-group invariant distributions for the eigenvalues. One 
then computes various associated distributions, e.g., eigenvalue spacing distribu- 
tions and their limits. For example, in the case of the Hermitian matrices the 
famous GUE-density appears as a limiting distribution in this way. 

At the very outset, one observes two major differences between the setting of 
disordered bosons and that of fermions. First, the necessity of dynamic stability 
imposes a strong condition on the operators being considered. Secondly, unlike 
the fermionic case where the model probability densities can be chosen to be 
invariant with repect to the full compact group of symmetries, in the bosonic 
case the group G = Sp2„(M) is noncompact and G-invariant densities are not 
available. 

In |LSZj an interesting class of bosonic emsembles {£, P) is introduced and the 
program which we sketched above for fermions is carried out. Here, as in |LSZj . 
ensemble densities which are invariant with respect to the (unique up to con- 
jugation) maximal compact subgroup K oi G are proposed. The usual choice 
is K = Un- Our hope here is that i^-invariant constructions which are closely 
related to the geometry and represention theory of the symplectic group G will 
also produce bosonic ensembles of quantum mechanical interest. We begin by 
explaining this approach for arbitrary n. Calculations are carried out in what 
might seem to be a toy model, i.e., the case of n = 1. Due to the fact that 
we are dealing with a group of Hermitian type, where the Poincare disk plays 
a fundamental role in understanding the associated symmetric space, this is an 
important special case in our future work. 



1. Generalities 

From the mathematical point of view there is a very interesting canonically de- 
fined neighborhood U{G) of the 0-section of the cotangent bundle T*M of the 
associated symmetric space M = G / K . In the terminology of ^FHW], to which 
we refer for background, it is called the universal domain associated to G. Its 
importance in representation theory, where it is most often called the crown of 
the symmetric space, has been underlined in numerous works. In the case at 
hand, where G is of Hermitian type, this is G-equivariantly and holomorphically 
isomorphic to the product B x B of the associated Hermitian bounded domain 
and its complex conjugate B ([BHHJ). It should be noted that B and B are 
biholomorphically equivalent as complex manifolds, but are not G-equivariantly 
biholomorphic. 
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Background on Hermitian symmetric spaces. — In order to emphasize 
the concrete nature of this situation let us provide a sketch of some basic in- 
formation about the Hermitian symmetric space B (see, e.g., |H] for details). 
A symmetric space is simply a (finite-dimensional) Riemannian manifold (M, g) 
with the additional property that at every p & M there exists an isometry cxp 
which is a symmetry at p in the sense that cXp = Mm and dap = — Idj^M- One 
shows that the connected component of the identity of the isometry group is a 
Lie group G which acts transitively on M and that the isotropy subgroup at a 
neutral point under consideration is a compact subgroup K. 

Based on curvature conditions and canonical decompositions of Lie groups, sym- 
metric spaces break up into products of symmetric spaces of three irreducible 
types. One such type is M = G/K where G is a simple (real) noncompact Lie 
group and i^' is a (unique up to conjugation) maximal compact subgroup. In 
some cases these symmetric spaces M have the structure of a complex manifold 
such that the group G acts as a group of holomorphic transformations. This 
situation can be group-theorectically characterized by the condition that K has 
a positive-dimensional center which in fact turns out to be In this case the 
metric can be chosen to be Hermitian and therefore one refers to M as a Her- 
mitian symmetric space and G being of Hermitian type. Note that the maximal 
compact subgroup K = Un oi G = Sp„(]R) satifies this condition. 

Noncompact Hermitian symmetric spaces of the type discussed above are nat- 
urally realized as distinguished open G-orbits in compact Hermitian symmetric 
spaces Z which are referred to as the compact duals of the noncompact Hermitian 
symmetric spaces. A simple example is the case where G = SU(n, 1) acts on the 
space Pn(C) of lines in C"^^. Here an associated symmetric space is the open 
G-orbit of negative lines. In general there are always exactly two G-invariant 
complex structures on a Hermitian symmetric space. In this case the "other" 
structure is space of positive n-planes in P„. 

Using a precise description of the action of the complex Lie group G^ on Z, 
one observes that the noncompact Hermitian symmetric space is contained as 
a bounded domain in a canonically determined (dense) open subset of Z which 
is biholomorphically equivalent to C™ where m := dimcM. Conversely, every 
(irreducible) bounded domain D in which is symmetric with respect to holo- 
morphic symmetries arises in this way. In our case we may view the Hermitian 
symmetric spaces B and B as open G-orbits in the space Z of Lagrangian sub- 
spaces of (C^"',Ci;), where u = dp A dq is the standard (holomorphic) symplectic 
structure. 

As a complex manifold X = B x B can also be G-equivariantly realized as a G- 
invariant open neighborhood of the (totally real) orbit G.zq of the neutral point 
in the affine symmetric space G'^/K'^. As mentioned above, it is G-equivariantly 
identifiable with a neighborhood of the 0-section in T*M. These identifications 
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allow US to consider natural G-invariant symplectic structures on X and the 
resulting moment maps /i : X — )■ g*. The most obvious of these are the restriction 
to X of the canonical symplectic structure on T*M and structures of the type 
uj := iddp, where p is a strictly plurisubharmonic G-exhaustion of X defined, 
e.g., by representation theory. 

One can show that with respect to any of these structures the generic G-orbits 
are coisotropic, i.e., that the /x- fibers are contained in the G-orbits. Thus an 
associated moment map transports a setting of geometric and representation 
importance to its image in an optimally controllable way. Our goal is 

to construct bosonic ensembles by equipping such an image with /i-direct image 
measures. However, for the symplectic structures mentioned above, p{X) is never 
contained in £. 

The fact that this happens, which was initially surprising to us, should be further 
investigated. Here, however, we move to a slightly different complex geometric 
situation by implementing the G-equivariant complex conjugation k, : B ^ B 
and thereby redefining X = B x B. Thus we still have the advantage of the 
symplectic geometry and representation theory of the original setup. In addition, 
in this new situation the moment map fi associated to any G-invariant Kahlerian 
structure on X satisfies fi{X) C S. Furthermore, /i is generically an open map 
and consequently its image is not pathological. The discussion of the /i-images, 
both before and after complex conjugation, is contained in ^ 

In the following section we compute fi for a particularly beautiful symplectic 
structure on the product X = A x A, where A is the unit disk in the complex 
plane. Here we use the identification of G with SU(1, 1). This is certainly a toy 
model. However, the bounded domain B contains a (flat) r-dimensional complex 
polydisk A^, defined by a choice of a maximal (abelian) subalgebra a in the p- 
part of a Cartan decomposition g = t ©p such that K.A^ = B. Therefore certain 
results can be immediately carried over to B by arguing one variable at a time 
in Ar- On the other hand, our goal of understanding asymptotic information 
concerning spectral invariants is of course not touched here. Nevertheless we 
do indicate how to compute the distribution function of the eigenvalue iu with 
oj > 0, and we compute this explicitly in the concrete case of the direct image of 
the uniform measure. 



2. The SU(1, l)-toy model 

Define the mixed signature Hermitian form (, )i i on by {z,w)i^i := ZiWi — 
Z2W2 and denote by G = SU(1, 1) its group of complex linear isometries. As usual 
g := Lie(G). The group has exactly three orbits on the projective space Pi(C): 
A := G.[0 : 1], A = G.[l : 0] and the real projective space Pi(M) = G.[l : 1]. 
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Note that the involution k : Pi — )■ Pi, [zi,Z2] [z2 : zi] defines a G-equivariant, 
antiholomorphic isomorphism k : A — > A. 

For moment-map considerations it is convenient to define the symmetric bihnear 
form 6 : X g ^ M by b{x, y) :— ^Tr{xy) which is invariant under the adjoint 
representation of G. Using the duahty defined by b we will regard the moment 
map as having image in g. We will use the following 6-orthogonal basis of g: 

If a moment map with values in g* is given by /i — /x^^* + /i^?7* + fi(C*: then the 
corresponding map with values in g is given by /i — —/i^C + A*??^ + A*cC- 

Kahlerian structure. — Here we make two observations which make a case 
for considering the B x B instead of S x S as the appropriate phase space. 

Proposition 2.1. — The image of a moment map /x : A x A — >■ g* which is 
defined either by the restricted coadjoint structure or by a Kdhler form uj = iddp 
of a G-invariant potential is contained in the complement of the cone 8 of elliptic 
elements. 

For a complementary statement for X = A x A we denote by £ the augmentation 
of either the positive or negative cone of elliptic elements by adding e g. 

Proposition 2.2. — The image of a moment map : A x A — > g* which is 
defined by a Kdhler form u = iddp of a G-invariant potential is contained in an 
augmented cone S with the diagonal being mapped to and the complement of the 
diagonal being mapped to the component of S which is contained in S. 

It should be noted that these observations can be translated to the higher dimen- 
sional setting by using strongly orthogonal roots and the polydisk slice mentioned 
above. 

Turning to the proofs of these remarks, we first note that in the case where the 
Kahler form is defined by an invariant potential, an associated moment map is 
given by 

(1) =: f^d^) = -Jiix){p) . 

Here J is the complex structure and ,^ denotes the vector field on the manifold 
which is associated to e g. Recall that since G is semisimple, the moment 
map is the unique equivariant map to g* with the Hamiltonian property dp^ = 
L^uj. Of course ^ is tangent to the level sets of the invariant function p and the 

moment map is measuring how p is growing along — J^. In this situation it is 

therefore important to determine the complex tangent spaces T^'^{p = p{x)} 
of the p-level sets. Assuming dp{x) ^ 0, this Cauchy-Riemann tangent space 
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is the 1-codimensional maximal complex subspace of the (real) tangent space 
Tx{p = p{x)}. 



Proof of Proposition \2.1\ Let 2:0 G A x i? be the base point so that M = G.zq = 
G/K is the unique symmetric space orbit in A x A. In the standard homogeneous 
coordinates this is defined by = 0. Choose a C p to be spanned by the 

matrix 77 above. Let S be the image of the half-open interval [0, |) by the map 
1 1—)- exp{tiri). Then S is an exact slice for the G-action on A x A (See |FHW J for 
this and other background used here.). This means that every G-orbit intersects S 
in exactly one point. The orbits Mp := G.p for p G S\{2o} are real hyp ersurf aces. 
The complex tangent spaces T^^Mp are calculated in a general setting in |FHWj 
(see p. 103). Applying this to the case at hand, we observe that these spaces are 
spanned by z>(p) and 6'z/(p) where u spans the +2 root-space defined by 77 and 6 
is the Cartan involution. In concrete terms = C,+C ^-^id = ^ — (. As we have 
seen above, it follows that p^ip) = Pevip) = 0. It follows that /i(S) is contained 
in the 77-axis where all isotropy groups are noncompact. In particular, the image 
of n is contained in the complement of S. 

It remains to prove the same result for the restriction of the canonical symplectic 
structure on T*M. Recall that in this case /i is defined by the Liouville form A; 
in particular, fi^{p) = ^{p){^{p))- However, if vr : T*M —)■ M is the canonical 
projection, then 7r^,(p)(,^) = and consequently fi({p) = for all p G S. Thus 
again we observe that Im(/i) has empty intersection with S. □ 



Proof of Proposition 2^ Here we use the slice S which is defined as the image of 
the half-open interval [0, 1) by the mapping t 1— ?■ {t, —t) in the affine coordinates 
{z,w) corresponding to ([1 : z], [1 : w]) in Z := Pi x Pi. The base point for these 
considerations is := (0,0). Just as above we compute the complex tangent 
space Tp^Mp of an orbit G.p = Mp for p G S \ {0}. For this observe that in 
this case p is in the open C^-orbit in Z. Therefore the complex tangent space 
TpZ can be represented as the quotient space Q^/Qp, where Qp is the (complex) 
Lie algebra of the G'^'-isotropy group at p. We note that this is conjugate to the 
subgroup of diagonal matrices in G = SL2(C). 

A direct calculation (See [S] p. 74) shows that 0^ consists of those matrices of 
the form 

(2) A=(l f ) with c G C. 

It is then immediate that t] and ( are linearly dependent modulo when regarded 
as elements of g^. Thus the real subspace of TpZ which is generated by 57 (p) and 
({p) is a complex line which by dimension arguments is the complex tangent 
space T^^Mp. 

Finally, by exactly the same argument as in the proof of the previous theorem it 
follows that Hr^ij)) = pdp) = for all p G S. Since /i = — /i^^ + + /^cC; 
follows that yuS = — /i^'C- Since is necessarily a critical point of p, it follows that 
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/i(0) = 0. Standard computations of rank(/x*) show that /^(p) 7^ for p e S \ {0} 
and the desired result follows. □ 
Actually we proved more than was stated in the above proposition. For future 
reference we state this here. 

Zusatz. The image //(S) is contained in the ^-axis {rj* — (^* — 0}. 



3. Poincare moment map 

Here we let dp : X ^ M.-^ be the distance function defined by the Poincare 
metric on A. We temporarily remove the diagonal from X where dp is no longer 
smooth, show that dp is strictly plurisubharmonic so that ou — idddp is a G- 
invariant Kahler form and compute the restriction to the slice S of the resulting 
moment map. As would be expected this map can be (continuously) extended to 
X with the diagonal being mapped to 0. 

Strict plurisubhcirmonicity of dp. — It is convenient to change coordinates 

on X, letting u = z + w and v = z — w where {z,w) arc the affine coordinates 
which are used above. Thus the complex disk L := {u = 0} contains the slice S 
as a radius. In particular L is transversal to the G-orbit of any point p G S \ {0}. 
Note that L is invariant by the diagonal 5'^-action defined by the element ^ G 5. 
It therefore follows that a function / on L is strictly plurisubharmonic if and only 
if f{e^) is strictly convex. Let us begin with this step which of course requires an 
explicit computation with dp. 



If 

ds{z,w) :-- 
is the Schwarz distance function, then 



z — w , 

1 — IVZ ' 



dp = log ]'^^f = tanh ^{ds) . 
L - ds 

Since tanh~^ is strictly convex, in order to prove the following it is enough to 
show that /(e^) is strictly convex where / is the restriction of ds to L. 

Lemma 3.1. — The restriction of dp to L is strongly suhharmonic. 

Proof. — Using the obvious coordinates. 



1 + ■ 

□ 

Now we turn to the "other" direction, namely the coordinate u. For this at each 
point of L we consider the complex curve 7 defined by 7(ti) = (t + ti, — t + ti). 
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Lemma 3.2. — Near m = the pullback 6{u) = dp{'y{u)) — dp{'y{0)) is nonneg- 
ative with 6{0) = and 6{u) > otherwise. Furthermore, 

(3) ^(0)>0. 

Proof. — As above it is enough to prove the positivity for dp replaced by ds- 
This is equivalent to 

\l-iu + t){u-t)\^ <{l + t^f 

for |m| sufficiently small and nonzero. This is done by elementary manipulations 
(See jS], p. 66-7). One can also prove ^ by making several estimates (Also see 
[S] , p. 67), but this can be proved in a more conceptual way: Since the curve 
7 touches the dp level set from outside, it is immediate that |^(0) > 0. If this 
derivative vanished at p, then the Levi-form restricted to the complex tangent 
bundle of the orbit G.p would vanish identically. It would then follow that G.p 
would be foliated by complex curves. This is for numerous reasons impossible, 
e.g., continuously moving the leaf through p away from p so that it no longer 
intersects the curve 7 yields a contradiction to Hurwitz's Theorem. □ 

Proposition 3.3. — The Poincare distance function is plurisubharmonic on X 
and is strictly plurisubharmonic outside of the diagonal. 

Proof. — The above two lemmas show that 

d^dp a^dp 

> and— ^ > 



at every point of L \ {0}. The fact that the curves 7 are tangent to the dp-level 
sets at every point of L \ {0} implies that vanishes along L and therefore 

the mixed derivatives vanish identically there as well. In other words the 
full complex Hessian is diagonalized along L \ {0} with positive entries along 
the diagonal. Since dp is only continuous along the diagonal of X, to show that 
it is plurisubharmonic there, we must prove that if /i : A — X is holomorphic 
with h{0) = 0, then dp o h is subharmonic. If the image of h is contained in the 
diagonal, then this is just the identically zero function. Otherwise we may assume 
that is the only point of the diagonal in the image and the result follows from 
the plurisubharmonicity of dp outside the diagonal and that, since dp{0) = 0, the 
mean value property is fulfilled at 0. □ 

Restricted moment map. — To keep the notation straight let Ex denote the 
slice S in X. The sign of the moment map chosen so that its image is contained 
in the augmented cone S, where S is the positive cone defined as G.He where 
'■= M^°^. By the above Zusatz we know that yu|Sx : Ex — )■ S^;. Since Ex and 
T,E are slices for the respective G-actions, for our purposes here it is enough to 
have an explicit description of the restricted moment map /u|Sx. This amounts 
to a description of the function fi^ along Ex- 
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Proposition 3.4- — If L is parameterized by z ^ {z, —z) and the ordered basis 
{iiViC} is used for Q, then 
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z 




1 - 


z 


2 



,0,0). 



Proof. — In this parameterization the restriction of the Poincare distance funtion 
to L is given by 

dp{z) = log 



By definition 



= -J^(^z)dp = - 



d_ 
ds 



dp{e''^z) = -j- 
s=o ds 



log 



1 + e 



-2s I 



-2s I 



and the desired result follows by a direct computation of the derivative. 



□ 



Corollary 3.5. — The moment map n : X ^ S is surjective. The diagonal in 
X is the fi-preimage of the vertex and the restriction of fi to the complement of 
the diagonal is a proper -bundle. Furthermore, fi has the coisotropic property 
that fi~^{fi{x)) C G.x for all x G X. 



4. Direct image of the uniform measure 

Let us conclude this note by explicitly computing the spectral density for the 
direct image measure P on E defined by the uniform measure u on X. It appears 
that other natural densities, e.g., those related to Gaussian type distributions 
will only be numerically computable. Closed form descriptions of the behaviour 
of these at the vertex of the cone should be possible. 
Note here that the matrices on the slice T^e are of the form 

with u > 0. The relevant distribution is therefore F{x) := P{u} < x). We will 
sketch the explicit computation of F here, referring to [S] for details. 
The mapping fj,\T,x '■ is given by x = = Let t = t{x) be the 

corresponding inverse value. Thus F{x) is the Euclidean volume of the region 
{dp < t{x)} = R{x). We compute this via fiber integration using the projection 
TT : X — 7- A on the first coordinate. For a £ A let A{a,x) be the Euclidean area 
of the fiber vr~^(a) fl R{x) =: D{a, x) so that up to a constant 

F(x) = / A(a,x)daAda. 

J A 

Let us begin by computing A{0,x). 

Lemma 4-1- — The fiber R{x) fl 7r^^(0) is a disk of radius u = u{t{x)) = 
in the TT-fiber {0} x A. 



10 



ALAN HUCKLEBERRY AND KATHRIN SCHAFFERT 



Proof. — Since the projection tt is G-equivariant, i?(a;) n7r~^(0) is a region in the 
TT-fiber which is bounded by an orbit of the G-isotropy group at in the 7r-image 

space. Since this isotropy group is acting on the fiber by rotations, this is a disk. 
Now the region R{x) fl L in L has the point pu = {u, —u) on its boundary. Thus 
if we apply the tranformation g{z,w) = {Tu{z),Tu{w)) with 

then Pu is mapped to the point (0, T„(— m)) which is on the boundary of the disk 
in question and the desired result follows. □ 

Since we are only dealing with the uniform measure on X, it follows that 



-^(^) = 77" / A{a,x)daAda. 

271" JaeA 



Now D{a, x) is the image of the disk D{0, x) by the transformation T^. Since 
preserves the Poincare metric, this is a hyperbolic disk, i.e., a disk with respect 
to dp. But hyperbolic disks are Euchdean disks. Thus, 

-^(^) = o / r"^ {a, x) da A da 

where r = r{a,x) is the Euclidean radius of D{a,x) which we now compute. 
By rotational invariance it is enough to compute r(a, a;) for a = s e [0, 1) on the 
positive real axis of A. Thus the transformation Tg stabihzes the interval / of real 
points of A. Since / is orthogonal to the boundary of D{0, x) and Tg is conformal, 
it follows that / is orthogonal to the boundary of D{s,x). Consequently, the 
length of / n D{s, x) is the diameter of D{s, x). This leads to the following fact. 

Lemma 4-2. — For s e [0, 1) the radius of D{s,x) is 

r{s,x) = -iUu) - U-u) = ^^-^ . 
Introducing polar coordinates and using rotational symmetry, one shows that 



sds 



for u — u{t{x)) as above. This integral can actually be computed in closed form: 

F(x) = 2—^ log(l -v?)^2-v? . 

Using the simple dilation x = | to clean up the numbers, one computes that 

;2 



.2 



X 



u 



1+^2 

Thus one has a rather simple closed form description of the distribution function 
F and the associated density /. 
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Proposition 4 '3. — In the variable x — ^ the eigenvalue distribution function 
F and its density f associated to the direct image of the uniform measure by the 
Poincare moment map /i : X ^ £ are given by 



p . ^ ' 1+^2 

and 

/ = — log 1 + x)- . 

x'^ x{i + x^y 

Numerical computation of / yields the following picture. 




Figure 1. Density function 



Near one has the power series representation 

k=2 

In particular f ^ x^ near 0. For x ~ oo the expression in the above proposition 
implies that / ~ ■ The expected value E(X) — xf{x)dx can actnally be 
explicitly computed as E{X) = |7r. However, as one observes from the estimate 
of / near oo the second moment E{X'^) is infinite. This phenomena does occur in 
nature, but here it is probably due to the fact that the original uniform measure 
on the phase space X is physically unreahstic. Gaussian densities, e.g., of the 
type e~^p would perhaps be of greater interest, but the resulting integrals are 
more complicated. 
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